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Corc-coUapse theory brings together many facets of high-energy and nuclear astrophysics 
and the numerical arts to present theorists with one of the most important, yet frus- 
trating, astronomical questions: "What is the mechanism of core-collapse supernova 
explosions?" A review of all the physics and the fifty-year history involved would soon 
bury the reader in minutiae that could easily obscure the essential elements of the phe- 
nomenon, as we understand it today. Moreover, much remains to be discovered and 
explained, and a complicated review of an unresolved subject in flux could grow stale 
fast. Therefore, in this paper I describe what I think are various important facts and 
perspectives that may have escaped the attention of those interested in this puzzle. Fur- 
thermore, I attempt to describe the modern theory's physical underpinnings and briefly 
summarize the current state of play. In the process, I identify a few myths (as I see 
them) that have crept into modern discourse. However, there is much more to do and 
humility in the face of this age-old challenge is clearly the most prudent stance as we 
seek its eventual resolution. 
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I. INTRODUCTION 

Stars are born, they live, and they die. Some, the 
most massive (> SM©), die explosively, spawning in the 
process neutron stars or "stellar-mass" black holes while 



littering the interstellar medium with many of the ele- 
ments of existence. But what is this process by which a 
star's multi-million year life is terminated abruptly and 
violently within seconds, then announced over months 
via the brilliant optical display that is a supernova explo- 
sion? Fifty years of theory, calculation, and observation 
have not definitively answered that question, though a 
vivid picture of the mechanism and terminal scenario of 
the dense core of a massive star is emerging. 

The solution to the puzzle of the mechanism of core- 
collapse supernova (CCSN) explosions involves more 
than just obtaining simulation explosions on supercom- 
puters. If this weren't the case, theorists would have 
solved this poser many times (and have!). Rather, the 
"solution" involves quantitatively explaining a host of 
astronomical facts that surround the supernova phe- 
nomenon. These include, but are not limited to: 1) 
the canonical explosion energy of ~10^^ ergs (defined as 
one "Bethe"), along with its putative distribution from 
^0.1 Bethes to ^10 Bethes. Explosion energy is poten- 
tially a function of progenitor mass, rotation rate, mag- 
netic fields, and metallicity. To date, no one has come 
close to achieving this central goal; 2) the residual neu- 
tron star mass and its distribution as a function of star. 
This involves more than simply noting that the Chan- 
drasekhar mass (Mch) of ^1.4 M0 is similar to the grav- 
itational masses of well-measured neutron stars (though 
this fact is relevant to zeroth-order) . The proto-neutron 
star (PNS) fattens by accretion during the respite before 
explosion, so the density and angular-momentum pro- 
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files in the progenitor core, the time of explosion, and 
the amount of fallback are all determining factors that 
are intimately tied to the mechanism and its unfolding. 
The branch map connecting progenitor to either neutron 
star or stellar-mass black hole final states is a related 
goal; 3) the nucleosynthetic yields as a function of stel- 
lar progenitor. Which stars yield how much of which 
elements is a combined function of a) the pre-explosion 
stellar evolution to the "onion-skin" structure of progres- 
sively heavier elements as one tunnels in to the central 
core and b) the explosion itself, which determines the 
mass cut and the degree of explosive nucleosynthetic re- 
processing. The yields of the elements between calcium 
and the iron peak (inclusive) are particularly sensitive to 
the explosion process, constituting as they do the inner 
cjccta; 4) the high average pulsar proper motion speeds. 
Radio pulsars are the fastest population of stars in the 
galaxy, with average speeds of ^350 km s~^, but rang- 
ing beyond ^1000 km . Asymmetries in the explosion 
itself and simple momentum recoil are natural culprits, 
but we don't have an explanation for their observed speed 
spectrum, nor for which progenitors give birth to the low 
proper-motion subclass of neutron stars bound as accret- 
ing X-ray sources to globular clusters; and 5) supernova 
explosion morphologies and ejecta element spatial distri- 
butions. Instabilities and asphericities in the explosion 
itself are compounded by instabilities during the propa- 
gation of the supernova shock wave through the progeni- 
tor star and the circumstellar medium to create a debris 
field that is anything but spherical and neatly nested. 
Even a qualitative identification of the signatures of the 
explosion process itself in the density and element distri- 
butions of the expanding supernova blast and subsequent 
supernova remnant (SNR) would be an advance. 

The literature on core-collapse supernova theory is 
vast and entails a fifty-year history of hydrodynamics 
and shock physics, radiative transfer, nuclear physics 
(at many junctures), neutrino physics, particle physics, 
statistical physics and thermodynamics, gravitational 
physics, and convection theory. The explosion signatures, 
along with those listed in the previous paragraph, are 
photon light curves and spectra, neutrino bursts, gravi- 
tational wave bursts, and meteoritic and solar-system iso- 
tope ratios. Thousands of researchers have at one time or 
another been engaged, many as careerists. To attempt 
to summarize or synthesize this literature, even with a 
focus on the theory of the mechanism, would be a gar- 
gantuan undertaking. Moreover, since the fundamental 
mechanism has not been satisfactorily and quantitatively 
demonstrated, such an ambitious review might seem pre- 
mature. 

Nevertheless, there have in the past been attempts to 
review core-collapse explosion theory, and some of these 
contain useful information and perspectives. All, how- 
ever, due to the inexorable evolution of the subject as 
researchers have struggled towards the ultimate goal of 



understanding, and due to the vastness of the task, are 
of limited scope and clearly trapped in time. This is 
natural. However, for those readers who desire a Cook's 
tour in the tradition of a standard, though helpful, re- 
view of the various aspects of core-collapse theory, I list 
here a sampling. Neither the samples, nor the sampling 
arc complete, and many of these papers were not writ- 
ten as reviews. For overall perspectives, I point to Bethe 
(1988;1990), Janka (2001,2012), Janka et al. (2007), Bur- 
rows et al. (1995,2007ab), Hcrant et al. (1994), Kotake, 
Sato, & Takahashi (2006), Kotake et al. (2012ab), Bur- 
rows (2000), and Mezzacappa (2005). For neutrino mi- 
crophysics, one can consult Dicus (1972), Bruenn (1985), 
Burrows, Reddy, & Thompson (2006), Tubbs & Schramm 
(1975), Freedman, Schramm, & Tubbs (1977), and Lan- 
ganke & Martinez-Pinedo (2003). For equation of state 
issues, there are Lamb et al. (1981), Lattimer (1981), 
and Lattimer & Swesty (1991). For massive star evo- 
lution and nucleosynthesis, good sources are Burbidge 
et al. (1957), Weaver, Zimmerman, & Woosley (1978), 
Woosley & Weaver (1986,1995), Thielemann, Nomoto, & 
Hashimoto (1996), Weaver & Woosley (1993), Nomoto et 
al. (1997), Woosley, Heger, & Weaver (2002), and Heger, 
Woosley, & Spruit (2005). For the connection with 
gamma-ray bursts, one can turn to Woosley & Bloom 
(2006) and for gravitational wave signatures there are Ott 
(2009), Kotake (2012), Miiller et al. (2004), and Miiller 
et al. (2012). There are many papers on the computa- 
tional issues specific to core collapse, but the papers by 
Arnett (1966,1967,1977), Imshennik & Nadczhin (1973), 
Bowers & Wilson (1982), Bruenn (1985), Maylc & Wilson 
(1988), Liebendorfer et al. (2001), Liebendorfer, Rampp, 
Janka, & Mezzacappa (2005), Burrows et al. (2000), 
Livne et al. (2004,2007), and Swesty & Myra (2009) 
arc collectively educational. Ph.D. theses by R. Mayle, 
A. Marek, and B. Miiller are particularly informative. 
Classics in the subject, though not particularly compre- 
hensive, include Bm-bidge et al. (1957), Colgate & John- 
son (1960), Hoyle & Fowler (1960,1964), Fowler & Hoyle 
(1964), Colgate & White (1966). Arnett (1966), LeBlanc 
& Wilson (1970), Wilson (1971), Mazurek (1974), Sato 
(1975), Bethe et al. (1979), and Bethe & Wilson (1985). 

To this set of "reviews," and in particiilar to the list of 
"classics," one could add many others. However, anyone 
well-versed in the papers on the above short list will be 
well-informed on most of supernova theory, if not com- 
pletely au courant. Nevertheless, there are many facets 
of theory, and various points of principle (some rather 
crucial) , that have been insufficiently emphasized and ar- 
ticulated in detailed research papers, which of necessity 
present a given, though sometimes narrow, result. The 
upshot has been some defocusing of the theory enter- 
prise and the accumulation of various myths that, while 
understood to be such by most practitioners, have at 
times confused the uninitiated. With this paper, I high- 
light an eclectic mix of important topics in core-collapse 
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supernova theory that I feel have not gotten sufficient 
"air time." In the process, some of the central features 
of modern supernova theory are identified. It is hoped 
that this collection of excursions, though idiosyncratic, 
will help sharpen collective understanding of the core is- 
sues to be tackled on the way to a complete and credible 
understanding of the explosion phenomenon. 



II. PHYSICAL CONTEXT OF CORE-COLLAPSE - 
BASIC SCENARIO 

A. Progenitors 

A central facet of a star's evolution is the steady de- 
crease with time of its core specific entropy. The loss 
by outward diffusion of the energy generated by nu- 
clear transmutation to progressively heavier elements (for 
which nucleons are more and more organized in the nu- 
cleus) naturally leads to these lower entropies. This may 
seem counter-intuitive, since core temperature increases 
during and between burning phases. However, core den- 
sity also increases, and this increase outpaces the tem- 
perature increase needed to maintain burning and, as ash 
becomes fuel, to ignite the next burning stages. In addi- 
tion, after the ignition of core carbon burning, the tem- 
peratures are sufficient to generate high fluxes of thermal 
neutrinos. These stream directly out of the core, accel- 
erating core evolution and, in the sense of entropy, its 
refrigeration. This phase also marks the evolutionary de- 
coupling of the stellar photosphere from the burning core, 
which then inaugurates its more rapid race to collapse. 

Entropy is relevant because it is a measure of random- 
ness. A low value bespeaks organization and electron 
degeneracy. Hence, stellar evolution leads to white dwarf 
cores. The more massive stars generate more massive 
cores, and those more massive than ~8 achieve the 
Chandrasekhar mass, at which point such cores are un- 
stable to the dynamical implosion (collapse) that inau- 
gurates the supernova. However, the "Chandrasekhar 
mass," and its core density profile when achieved, depend 
upon the electron fraction (Yg) and entropy profiles. The 
effective Chandrasekhar mass is not its canonical zero- 
entropy, uniform- Ye value of 1.456(2Ye)^ M©. Entropy 
and Ye profiles are functions of the specific evolutionary 
paths to instability, in particular the character of con- 
vective shell burning and the ^^C(a,7)^^0 rate. More 
massive progenitors evolve more quickly, and, therefore, 
don't deentropize as much by neutrino losses before col- 
lapse. They have higher entropies, resulting in additional 
pressure beyond that associated with zero-temperature 
electron degeneracy to support more mass. Therefore, 
the core mass necessary to go unstable is increased. Note 
that electron capture on nuclei alters Ye, and the rates 
of electron capture depend sensitively upon isotope and 
density. 



The result is a "Chandrasekhar mass" that could vary 
from ~1.25 to ~2.0 M©, depending upon progenitor mass 
and evolutionary details (Figure [TJ . The latter are not 
necessarily perfectly captured by current models. How- 
ever, the trend seems roughly to be that more mas- 
sive progenitors have 1) larger effective Chandrasekhar 
masses and 2) envelopes in which the mass density de- 
creases more slowly and, hence, that position more mass 
around the core. Figure [2] shows the core density profiles 
for various theoretical initial progenitor masses. While 
most massive-star supernova progenitor cores evolve to 
iron peak elements and a Chandrasekhar-mass "iron 
core" at the center of the canonical "onion-skin" struc- 
ture of progressively lighter elements from the inside 
out, the least massive progenitors (perhaps with ZAM^ 
masses of ~8. 0—9.0 M©) are thought to end up as "O- 
Ne-Mg" cores (Nomoto & Hashimoto 1988). Such cores 
might have very tenuous outer envelopes which, given 
current thinking, might result in underenergetic (^10^*^ 
ergs = 0.1 Bethe) neutrino- wind-driven explosions (Bur- 
rows 1987; Kitaura et al., 2006; Burrows et al. 2007c). 
In both circumstances, the ashes from shell burning are 
responsible for fattening the inner core to the effective 
Chandrasekhar mass. Importantly, the more massive 
progenitors have slightly lower central densities and tem- 
peratures at collapse. The higher densities of "0-Ne- 
Mg" cores result in higher electron capture rates, but the 
higher entropies of the more massive progenitors result in 
a greater softening of the EOS due to photodissociatiorj^ 
Both processes facilitate the achievement of the Chan- 
drasekhar instability, which once achieved is dynamicaj^ 
Therefore, in a real sense how the Chandrasekhar insta- 
bility is achieved is secondary, and collapse proceeds sim- 
ilarly in both cases. Some would distinguish "electron- 
capture supernovae" as a different species of supernova. 
However, this is really not justified. The true difference 
is in the envelope density profile (Figure |2|. The lowest- 
mass massive progenitors have very steep density pro- 
files outside the inner core. These translate into sharply 
dropping accretion rates onto the proto-neutron star after 
bounce and before explosion. The dynamical differences 
and outcomes of "iron-core" or "0-Ne-Mg" collapse are 
more dependent, therefore, on the different outer core 
mass and density profiles, where the differences are ex- 
pected to translate into real differences in explosion en- 
ergy, kick speeds, residual neutron star masses, and opti- 
cal displays. Since neither core type contains much ther- 
monuclear fuel, unlike in the Type la supernova case of 



Zero- Age-Main-Sequence 

which decreases the effective 7 of the gas below the critical value 
of 4/3 

The collapse time to nuclear densities is then no more than ~350 
milliseconds, whatever the progenitor mass. An approximate 
characteristic dynamical time is very roughly ^40 ms/y/pio, 
where pio is in units of 10^*^ gm cm~^. 
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critical carbon-oxygen white dwarfs, burning does not in- 
hibit collapse to nuclear densities ('--^2.6x 10^'* gm cm~'^). 
However, the reader should note that, due to carbon and 
neon shell flashes, stellar evolution simulations up to the 
edge of collapse are very difficult for the lowest-mass mas- 
sive stars and only one group has provided a massive-star 
model below 9.0 Mq (Nomoto & Hashimoto 1988). 

I reiterate that the initial core profiles, along with 
whatever initial rotation may be present in the core at 
collapse, must determine the spread in outcomes of col- 
lapse. In a very real sense, "progenitor is destiny," a 
mapping complicated only by the randomness associ- 
ated with chaotic turbulence and instability dynamics 
and with the unknown initial perturbation spectrum im- 
posed by pre-collapse convective burning (Meakin & Ar- 
nett 2006,2007ab). Due to such stochasticity, it is ex- 
pected that a given progenitor star and structure will 
give rise to a distribution of outcomes (energies, proto- 
neutron star masses, kick speeds and directions, etc.), 
with a "a" that remains very much to be determined but 
that many hope (without yet much justification) will be 
small. The task before theorists is to determine the pro- 
genitor/supernova mapping. On the observational side, 
Smartt (2009) has recently attempted such a mapping 
by identifying progenitors in archival data to a handful 
of core-collapse supernovae. His preliminary finding that 
no progenitor to a supernova with such archival data is 
more massive than ^16 Mq is intriguing, but will need 
further investigation to confirm or refute. 

With this background, core collapse proceeds (in 
theory) similarly for all progenitors. The inner 
^5000—10,000 kilometers is the most relevant. Its dy- 
namical time is less than a second, while that of the rest 
of the star (with a radius of ~10^— 10^ km), most of which 
is comprised of hydrogen and helium, is hours to a day. 
Therefore, the dynamical inner core is decoupled from 
the outer shells and it is only when the supernova shock 
wave reaches them that they too participate in the explo- 
sive dynamics, but then only as shocked spectators. The 
initial core densities and temperatures are ~6— 15x10^ 
gm cm~'^ and ~6— 10x10^ Kelvin, respectively. The ini- 
tial core entropies are ^0.7-1.2 fc^/baryon, where ks is 
Boltzmann's constant, while the entropies in the outer 
fossil convective silicon and oxygen burning shells jump 
to ^2—5. There is a corresponding abrupt decrease in 
density at these shell boundaries, as well as increases in 
Ye, whose initial core values are ~0.4— 0.43. The outer 
shells are at lower densities, for which the electron cap- 
ture rates are low and Ye is very near 0.5 (Figurc[T]). The 
pressures are dominated by degenerate relativistic elec- 
trons, with slight thermal and Coulomb corrections, and 
for most progenitors most of the baryons are in nuclei 
near the iron peak. If we are dealing with an "0-Ne- 
Mg" core, burning on infall will rapidly convert it into 
an "iron-peak" core. For an iron core, the nuclei are 
in nuclear statistical equilibrium, which is a Saha equi- 



librium predominantly between nucleons in and out of 
nuclei, alpha particles, and the nuclei. Such a "chemi- 
cal" equilibrium is described, for a given nuclear model, 
by temperature (T), mass density (p), and Ye. 



B. Collapse 

The instability that is collapse occurs because the aver- 
age adiabatic 7 in the core is at or below 4/3. Photodis- 
sociation and domination by relativistic electrons guar- 
antee this. As collapse proceeds, T and p both rise. With 
the temperature increase, more nucleons evaporate from 
nuclei. If the number of free nucleons were to rise signifi- 
cantly, since they are non-relativistic ideal gases with a 7 
of 5/3, collapse would be halted and reversed before nu- 
clear densities were achieved. This was the supposition in 
the early 1970s. However, as Bethe et al. (1979) and oth- 
ers have shown, the increase in temperature of the nuclei 
populates excited nuclear states, which represent many 
degrees of freedom. The upshot is a significant increase 
in the specific heat and the regulation of the temperature 
increase, since energy that would otherwise be channeled 
into kinetic degrees of freedom is redirected in part into 
these excited states. The result is not only moderation 
in the temperature increase during collapse, but also in 
the production of non-relativistic free nucleons, thereby 
maintaining the domination of the pressure by the rel- 
ativistic electrons and preserving the nuclei during col- 
lapse. The consequence is collapse all the way to nuclear 
densities, at which point nuclei phase transform into nu- 
cleons that at such densities experience strong nuclear 
repulsion, severely inhibiting further compression. 



C. Bounce 

Within less than a millisecond, this stiffening of the 
equation of state halts and reverses collapse. During col- 
lapse, since the central speed and outer core speeds must 
be zero and low (respectively) the peak collapse speed is 
achieved in the middle of the collapsing core. Early dur- 
ing collapse, this results in a separation into a subsonic 
inner core, which collapses almost homologously (v oc r) 
and as a unit, and a supersonically infalling outer core. 
The peak speed of the inner shells of the outer core are 
roughly a constant large fraction of free-fall. Therefore, 
when the inner core achieves nuclear densities and re- 
bounds, because these two regimes are out of sonic con- 
tact, the subsonic inner core bounces as a unit and as 
a spherical piston into the outer core, which is still col- 
lapsing inward, thereby generating a shock wave at the 
interface (Figure |3|. This is the supernova shock in its 
infancy. 

During collapse, increasing density and temperature 
result in increasing electron capture rates on both nu- 
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clei and free protons, with the resuhing decrease in Yg. 
Current thinking is that capture on nuclei predominates 
(Langanke & Martinez-Pinedo 2003), but since the "iron 
peak" shifts with the increase in p and the decrease in Y,, 
to higher and higher atomic weights, exotic isotopes, for 
which we have no data, quickly dominate. However, this 
uncertainty alters the progress of collapse only slightly, 
since gravitational free fall bounds collapse speeds, what- 
ever the capture rates. The result is that pre-bounce 
collapse is universal in character, requiring (after peak 
speeds have achieved ~1000 km s~^) between ~150 and 
~350 milliseconds to achieve bounce. More importantly, 
rewinding collapse from bounce yields almost universal 
trajectories relative to bounce time. 



D. Trapping 

As indicated, electron capture during collapse lowers 
Ye, but it also produces electron- type neutrinos (i/g) at 
progressively increasing rates. The average energy of 
these j/gS increases with density and temperature. Since 
the cross section for scattering off nuclei by the coherent 
process (Preedman 1974) increases roughly as the square 
of neutrino energy and the densities are soaring at a 
rapid rate, the mean-free-paths for i/g-matter interactions 
are fast decreasing. When these mean-free-paths become 
comparable to the scale-height of the matter density or 
when the outward diffusion speed of the V(,s equals the 
infall speed, then the z/gS are trapped in the flow (Mazurek 
1974; Sato 1975). After trapping, electron capture is bal- 
anced by Vf. capture to establish chemical equilibrium at 
a given electron lepton fraction. 

Since this condition is achieved not long after the cen- 
tral density is ~10^^ gm cm~"^, trapping of electron lep- 
ton number and z/gS happens before much electron cap- 
ture has occurred and has a profound c^ffcct on collapse 
and subsequent evolution. Trapping locks electron lepton 
number, electrons, and i^gS in the core for many seconds, 
depending on the mass shell. The trapped z^gS are com- 
pressed significantly, but at low entropy and conserving 
lepton fraction. The latter settles near 0.30 (partitioned 
between electrons [^0.25] and j^gS [^0.05]), instead of 
0.03 (all electrons), which it might have been without 
trapping. The result is the channeling by compres- 
sion of a significant fraction of the gravitational energy 
of collapse into degenerate i^gS and electrons, with cen- 
tral chemical potentials and Fermi energies of ~100— 300 
MeV. Hence, trapping is not merely the increase of the 
Vf, optical depth far beyond values of one, but the preser- 
vation at high values (not far from the initial value of 
^0.43) of the lepton fraction and electron numbers in the 
core and the production of a sea of degenerate v^, whose 
average energy, by dint of compression after trapping, is 
high. If it weren't for lepton trapping, the v^s would be 
thermal with average energies at bounce of ~30 MeV, 



optical depths of hundreds, and energy diffusion times of 
~50-100 milliseconds. Instead, their inner core average 
energies are 100-300 MeV, their optical depth to infinity 
from the center is >10^, and the energy diffusion time 
out of the core is many seconds. The latter has been 
boosted by the compressional increase in the average fg 
energy, for which the z/g-matter cross sections are much 
larger than they would otherwise be. 

By keeping Yg high, trapping thwarts the rapid evap- 
oration of neutrons from the nuclei that would otherwise 
be too neutron-rich to retain them. Therefore, both the 
increased specific heat due to excited nuclear states and 
trapping ensure the survival of the nuclei and the main- 
tenance of the effective 7 below 4/3 until nuclear densi- 
ties are achieved. If it weren't for trapping and excited 
nuclear states, bounce would occur at sub-nuclear densi- 
ties, the average core electron neutrino energy would be 
much lower, and fg-matter mean-free-paths at bounce 
would be mTich longer. The lower bounce densities and 
longer mean-free-paths would have translated directly 
into much shorter deleptonization and cooling times for 
a proto-neutron star. Hence, trapping is centrally im- 
portant in explaining the long duration of the neutrino 
burst from SN1987A (Kamioka II - Hirata et al. 1987; 
1MB — Bionta et al. 1987). Before neutrino trap- 
ping was recognized, the energy diffusion time out of 
the proto-neutron star was thought to be ^100 millisec- 
onds, a factor of ^100 shorter than observed. Moreover, 
since the total reservoir of energy radiated is fixed at 
the gravitational binding energy of a cold neutron star 
(^gm! _ 3 X 10''^'''ergs = 300 Bethes), a short dura- 
tion would have implied a higher average neutrino en- 
ergy (>50 MeV) than measured by Kamioka (~15 MeV 
— for the S'es). Therefore, and importantly, the lower 
measured energy and higher measured duration are di- 
rect consequences of electron neutrino (lepton number) 
trapping, a fact not widely appreciated. 

Trapping halts the short-term decrease in Yg during 
collapse, but, as mentioned, there is still a slight de- 
crease to ^0.25—0.35. The magnitude of this decrease 
is enhanced by inelastic scattering of i/gS off electrons 
(Mezzacappa & Bruenn 1993abc). Like Compton scat- 
tering, the capture-produced v^a are downscattered in 
energy and at lower energies their Freedman scattering 
cross sections are lower. The result is a slightly slower in- 
crease in the optical depth during infall, and this results 
in slightly delayed trapping at lower YgS. The magnitude 
of the effect is ~0.02-0.03. 

The trapped Yg sets the mass scale of the 
homologously-coUapsing inner core at approximately the 
associated Chandrasekhar mass (Yahil 1983; Burrows & 
Lattimer 1983). Since \ich is proportional to Yg^, this is 
~0.5— 0.6 M0. Therefore, the bounce shock first emerges 
at the sonic point near -^0.5— 0.6 Mq (~10-20 km) be- 
tween the inner core and the outer supersonic mantle in 
an optically thick region. It then propagates down the 
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density gradient, entropizing the matter it encounters, 
dissociating the nuclei into nucleons, and reaching lower 
Vf. optical depths within ~1 millisecond. At this point, 
the copious sea of z^eS, newly-liberated by electron cap- 
ture via the super-allowed charged-current capture pro- 
cess on newly-liberated protons, "breaks out" in a burst 
that lasts ~10 milliseconds (Figure |4]). The luminosity of 
this electron neutrino breakout burst is within an order- 
of-magnitude of the total optical output of the observable 
Universe (~3 — 4 x 10^^ ergs s~^), and is the first, most 
distinctive, feature of the supernova neutrino emission 
process. It is a firm prediction of generic core-collapse 
supernova theory, and if it doesn't exist, then much of 
the supernova theory developed in the last ^-^35 years is 
wrong. 



E. The Problem 

The direct mechanism of explosion posits that this 
bounce shock continues unabated outward into the star 
and is the supernova. However, both simple theory and 
detailed numerical simulations universally indicate that 
the burst and photodissociation of the infalling nu- 
clei debilitate the shock wave into accretion within ~5 
milliseconds of bounce. In a very real (though approxi- 
mate) sense, the subsonic inner core and shocked man- 
tle together execute a coherent harmonic oscillation that 
is near-critically damped. The shock acts like a black- 
body absorber of sound (the bounce pressure pulse), 
and the breakout neutrinos and photodissocciation do 
the rest. The breakout neutrino burst directly saps the 
shock of energy, while photodissociation by the shock 
redistributes shock energy from the kinetic component 
(and, hence, pressure) to pay the nuclear binding energy 
penalty. Contrary to common lore, photodissociation is 
not a loss of energy — the energy is still there and could 
be recovered with recombination. Rather, photodisso- 
ciation softens the equation of state by lowering the 7 
and raising the specific heat. The result is a less efficient 
conversion of infall kinetic energy (whose original source 
is gravity) into pressure. One should note that if elec- 
tron capture and neutrino transport are both artifically 
turned off during and after collapse, but a realistic EOS 
with photodissociation is employed, the direct mecha- 
nism can be shown to work for many of the published 
progenitor models. However, even then, with Yg frozen 
at its initial value and no neutrino burst or losses, the 
energy of the explosion is never higher than a few tenths 
of a Bethe, not accounting for the need to overcome the 
gravitational binding energy of the rest of the star above 
a "canonical" neutron star mass cut of ^1.5 (Figure 
[5]). Something more is needed. 

Therefore, since circa 1980 theorists have been pre- 
sented with a stalled accretion shock at a radius near 
~100-200 km and have been trying to revive it. This was 



and is an unsatifactory state of affairs. Supernova rates, 
O- and B-star death rates, and neutron star birth rates 
all suggest that most massive stars explode as supernovae 
and leave neutron stars. The fraction that branch into 
the stellar-mass black hole channel is unknown. The frac- 
tion of times stellar-mass black hole formation is accom- 
panied by a supernova is unknown. If the shock is not 
revived and continues to accrete, all cores will collapse 
to black holes. The serious revival mechanisms include 
the delayed neutrino mechanism (Wilson 1985), magne- 
tohydrodynamic (MHD) bipolar jet production (LeBlanc 
& Wilson 1970; Bisnovatyi-Kogan et al. 1976; Symbal- 
isty 1984; Burrows et al. 2007d; Takiwaki & Kotake 
2011) (requiring very rapid rotation), and the acoustic 
mechanism (Burrows et al. 2006,2007b). The spherical 
delayed neutrino mechanism works only weakly, and then 
only as a wind (Burrows 1987; Kitaura et al. 2006; Bur- 
rows et al. 2007c) for one non-representative progenitor 
(Nomoto & Hashimoto 1988), but the multi-dimensional 
variant is considered the frontrunner in the generic pro- 
genitor case. However, using the most advanced multi-D 
numerical codes and incorporating detailed microphysics, 
no one has yet been able to explain or reproduce any of 
the five "facts" of core-collapse supernova listed at the 
beginning of the introductiorj^ 

The stalled shock continues to dissociate infalling nu- 
clei into nucleons. Electron capture onto the newly- 
produced nucleons at densities of ~10^ — 10^° gm cm~^ 
and with temperatures of ^1.5 — 2.5 MeV rapidly lowers 
the post-shock Yg to values of ^0.1-0.2, creating the char- 
acteristic Ye trough seen in all detailed simulations (Fig- 
ure [6]) . Due both to electron capture behind the shock 
and neutrino cooling from the neutrinospheref]^ neutrino 
and electron-lepton-number losses continue apace. The 
proto-neutron star becomes more and more bound. 

Therefore, the object of current intense investigation 
is a quasi-static proto-neutron star (Burrows & Lattimer 
1986; Figure[7| with a baryon mass of ~1.2— 1.5 (de- 
pending upon progenitor), an initial central density at or 
above ~4.0 x lO^'' gm cm^"^ and rising, a central temper- 
ature of ~10 MeV (Figure [s]) and rising, an inner core 
entropy of ~1, and a shocked outer mantle entropy with 
a peak value ranging near ^-^6— 15 (Figure|9|, all bounded 
by an accretion shock wave stalled near ^100—200 km. 
During the delay to explosion, the accretion rate (Af) 
decreases from ~1.0 to 0.3 M0 per second (for the most 
massive progenitors) to from ~0.4 to 0.05 per sec- 
ond (for the least massive progenitors with iron cores). 
The associated accretion ram pressure, post-shock elec- 



^ except, perhaps, pulsar kicks: Scheck et al. 2004, 2006; Wong- 

wathanarat et al. 2010; Nordhaus et al. 2011 
^ at average radii of ~20-60 km, densities of ~10^^ to ~10^^ gm 

cm~^, and temperatures of ~4— 7 MeV of all six neutrino types 

(ve, Pe, Vii, i^fi, v-r, and Pr) 
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tron capture, and deep gravitational potential well are 
major impediments to explosion. If there were no accre- 
tion, neutrino heating would explode the proto-neutron 
star mantle immediately. The rate of accretion and its 
evolution are functions strictly of the progenitor mass 
density profile just prior to collapse, which is a central 
determinant of the outcome. 



III. THE CURRENT STATUS OF CORE-COLLAPSE 
SIMULATIONS 

There has been palpable progress in the development 
of techniques and tools to address the core-collapse prob- 
lem in the last thirty years, but the current status of 
the theory for the mechanism and the systematics of 
core-collapse explosions is ambiguous, if not confusing. 
Wilson (1985), in a pioneering paper and using a spher- 
ical code, obtained a neutrino-driven explosion after a 
short post-bounce delay, but only if he included a mixing- 
length algorithm to mimic doubly-diffusive "neutron- 
finger" instabilities that dredged up heat and thereby 
enhanced the driving neutrino luminosity and heating in 
the "gain region" (Bethe & Wilson 1985). The gain re- 
gion is the region behind the shock in which there is net 
neutrino heating. Without this boost, he did not ob- 
tain explosions. However, Bruenn & Dineva (1996) and 
Dessart et al. (2006) have shown that such instabili- 
ties don't arise in proto-neutron stars. More recently, 
Kitaura et al. (2006) and Burrows et al. (2007c) ob- 
tained weak ID (spherical) explosions (^10^" ergs) via 
a neutrino-driven wind for the lowest mass progenitor 
in the literature (~8.8 Mq). These are the only mod- 
els credibly shown to explode in ID (Mezzacappa et al. 
2001). This low explosion energy comports with the in- 
ference that the energy of explosion may be an increasing 
function of progenitor ZAMS mass (Utrobin & Chugai 
2009). Moreover, explosion occurred before convective 
overturn instability could be obvious. This result's ID 
simplicity lends credence to these state-of-the-art simu- 
lations using the best physics. However, this progeni- 
tor has a uniquely steep density gradient just exterior 
to the core (and, hence, almost no accretion tamp), and 
such progenitor density profiles and ID explosion behav- 
ior are seen in no other circumstance. This highlights 
an important point. There is not just one "core-collapse 
supernova problem," but many. The properties of presu- 
pernova stars vary and so too will the explosion model. 
More massive stars are more difficult to explode and at 
some point probably make black holes (Ugliano et al. 
2012; O'Connor & Ott 2011; Fryer et al. 1999). 

Shock-imposed negative entropy gradients at bounce, 
neutrino heating from below, the standing-accretion- 
shock- instability ("SASI"; Blondin, DeMarino, & Mez- 
zacappa 2003; Foglizzo, Scheck, & Janka 2006), and neg- 
ative lepton gradients in the inner core all render the 



shock-bounded environment of the proto-neutron star 
unstable and turbulent, severely breaking spherical sym- 
metry in the general case. Indeed, it was shown some 
time ago that multi-dimensional instabilities obtain and 
are probably central to the core-collapse supernova mech- 
anism in most cases (Burrows & Fryxell 1992; Herant et 
al. 1992,1994; Burrows, Hayes & FryxeU 1995; Janka & 
Miiller 1996; Blondin, DeMarino, & Mezzacappa 2003; 
Fryer & Warren 2002,2004; Foglizzo, Scheck, & Janka 
2006). The salient processes in the neutrino heating 
model may be the enhancement of the neutrino-matter 
heating efficiency, the turbulent pressure, and the en- 
largement of the gain region, among others. Neutrino 
heating in the shocked mantle exterior to the neutri- 
nospheres is probably the foremost driver of convec- 
tion (Burrows, Dolence, & Murphy 2012; Murphy et al. 
2012). This is an unexceptional conclusion, since in the 
neutrino heating mechanism, it is the driver of explosion 
and the major source of explosion energy. 

Using detailed neutrino transport and 2D hydro, 
Marek & Janka (2009) observed an explosion of a 15- 
Mq progenitor via the turbulent neutrino heating mech- 
anism. However, it seems underpowered (<10^° ergs) for 
this "average" progenitor, and explosion was not seen for 
their higher-resolution run. Recently, Miiller (2011) and 
Miiller, Janka, & Marek (2012) obtained explosions in 2D 
of both 11.2 Mq and 15-M0 progenitors using state-of- 
the-art neutrino transport with conformal general rela- 
tivity (emphasizing the importance of the latter; see also 
Kuroda, Kotake, & Takiwaki 2012), but the explosion 
energies were ~2.5xl0^^ and ^lO^^ergs, respectively, at 
the end of the calculations. When account is taken of the 
binding energy of the overlying star (Figure [5]) these en- 
ergies could be negative, potentially undermining these 
state-of-the-art 2D simulations as viable supernova ex- 
plosion models. However, the asymptotic energies re- 
main to be determined and it should be acknowledged 
that these are the best 2D simulations of collapse yet per- 
formed. However, as in Marek & Janka (2009), though 
the hydro was 2D, the transport in the Miiller (2011) and 
Miiller, Janka, & Marek (2012) papers was done multiple 
times in ID along numerous radial rays. The so-called 
"ray-by-ray" formalism uses the temperature (T), den- 
sity (p), and electron fraction (Ye) profiles along a given 
radial ray to feed a ID spherical transport calculation, 
that then provides local heating rates. Since vigorous 
convection in 2D leads to significant angular variation in 
T, p, and Yg, this approach exaggerates the variation 
in the neutrino heating with angle. Ott et al. (2008), 
performing the only multi-group, multi-angle neutrino 
transport calculations ever done in 2D, have shown that 
the integral nature of transport smoothes out the radi- 
ation field much more than the matter fields. This sug- 
gests that local neutrino heating rates obtained using the 
ray-by-ray approach might unphysically correlate with 
low-order mode shock and matter motions, and could be 
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pumping them unphysically. This may contribute to the 
perception that the £ = 1 dipolar shock osciUation seen in 
such 2D simulations is essential to the turbulent, multi-D 
neutrino-driven mechanism. 

Burrows et al. (2006,2007b) performed 2D radia- 
tion/hydro multi-group, flux-limited simulations and did 
not obtain neutrino-driven explosions. Their algorithm 
was Newtonian and did not include velocity-dependent, 
general-relativistic, or inelastic scattering terms in the 
transport. Nevertheless, ~1 second after bounce they ob- 
served vigorous inner core g-mode oscillations that gener- 
ated a sound pressure field sufficient to explode the enve- 
lope. However, the energy of explosion via this "acoustic" 
mechanism was very slow to accumulate, reaching only 
^^10^° ergs after more than ^^0.5 seconds. In addition, 
Weinberg & Quataert (2008) suggest that the amplitude 
of these g-mode oscillations could be severely diminished 
by a non-linear parametric resonance that could bleed en- 
ergy into very short-wavelength modes, dissipating this 
g-mode oscillation thermally. Such daughter modes are 
too small to simulate with current grids, but the weak- 
ness of the consequent explosions and the possibility of 
an important additional damping mechanism make this 
"acoustic" solution sub-optimal. Moreover, if the neu- 
trino mechanism in any form obtains, it would naturally 
abort the acoustic mechanism. 

Rapid rotation with magnetic fields should naturally 
lead to vigorous explosions and this MHD mechanism has 
a long history (LeBlanc & Wilson 1970; Symbalisty 1984; 
Akiyama et al. 2003; Akiyama & Wheeler 2005). The 
free energy of differential rotation available at bounce 
if the progenitor core is rapidly rotating {Pi ~ 1 — 2 
seconds]]^ is a potent resource, naturally channeled into 
bipolar jet-like explosions (Burrows et al. 2007d). Fig- 
ure [TO] portrays such a magnetic explosion. A variant of 
this mechanism has been suggested for gamma-ray bursts 
(MacFadyen & Woosley 1999). However, pulsar spin data 
indicate that the generic progenitor core spin rates must 
be rather low, and that at most ^1% of collapses can 
be via such an MHD mechanism. Hence, MHD-driven 
explosions can't be the generic core-collapse supernova 
channel. However, the neutrino heating mechanism may 
not be able to provide explosion energies above a few 
Bethes, and this suggests that the best explanation for 
the rare energetic "hypernovae" might be MHD power 
due to very rapidly rotating cores. 

There are numerous ambiguities and problems with the 
current generation of 2D simulations that are not repro- 
ducing the core-collapse supernova phenomenon. One 
symptom (or feature) of the problem in 2D may be the 
long time to explosion seen by, e.g., Marek & Janka 



translating into a core rotating with ~1.5- to 4-millisecond peri- 
ods 



(2009) , Suwa, et al. (2010), Miiller, Janka, & Marek 
(2012), and Burrows et al. (2006). Waiting -400-1000 
milliseconds to explosion may, certainly in the case of 
the neutrino-heating mechanism, "waste" the neutrinos 
emitted in the interval after bounce. However, Nature is 
not 2D and it could be that an impediment to progress 
in supernova theory over the last few decades is a lack 
of access to codes, computers, and resources with which 
to properly simulate the collapse phenomenon in 3D. 
This could explain the agonizingly slow march since the 
1960's towards demonstrating a robust mechanism of ex- 
plosion. The difference in the character of 3D turbulence, 
with its extra degree of freedom and energy cascade to 
smaller turbulent scales than are found in 2D, might re- 
lax the critical condition for explosion (Burrows & Goshy 
1993). Indeed, Murphy & Burrows (2008) have shown 
that the critical condition is more easily met in 2D than 
in ID and, recently, Nordhaus, Burrows, Almgren, & Bell 

(2010) conducted a similar parameter study comparing 
ID, 2D, and 3D and found that going from 2D to 3D 
could lower the threshold for explosion still more. How- 
ever, Hanke et al. (2012) have called this conclusion into 
question. Nevertheless, these developments suggest that 
it is a few tens of percent easier to explode in 3D than in 
ID, and that full 3D simulations, but with competitive 
multi-group neutrino transport, might be needed to prop- 
erly address this long-standing problem in computational 
nuclear astrophysics. Such "heroic" 3D simulations will 
be very computationally challenging, but arc the future 
(Takiwaki, Kotake, & Suwa 2012). 

IV. IMPORTANT FEATURES OF SUPERNOVA THEORY 

Following my description of the core-collapse scenario 
and my brief summary of the current status of the numer- 
ical theory, I now embark upon discussions of select top- 
ics, that though important, are often ignored, assumed, 
or misrepresented. However, I attempt merely to provide 
simple, yet useful, insights into basic supernova physics 
— rigor is not here my goal. In the process, I highlight 
some of the central themes and myths of core-collapse 
supernova theory. Since I conduct this survey in lieu of 
a final resolution of the supernova problem, the reader 
is encouraged to retain an open mind, and forgiven for 
retaining a critical one. 



A. General Themes 

Frequently missing in general discussions of core- 
collapse supernovae is that they are gravitationally- 
powered. Nuclear burning during explosive nucleosyn- 
thesis of the outer mantle after the explosion is well along 
might contribute at most —10% of the blast energy. A full 
solar mass of oxygen and/or carbon would have to burn 
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to iron peak to yield one Bethc. Given all extant progen- 
itor model profiles, much less than that amount of fuel is 
close enough to the core to achieve by shock heating suf- 
ficient temperatures (>4 x 10^ K). Moreover, before ex- 
plosion, any infalling fuel will be burned uselessly during 
collapse, and the ashes will be dissociated by the stalled 
shock and then biiricd in the core. The neutrino energy 
emitted from the core and absorbed in the proto-neutron 
star mantle that is required for the neutrino-powered 
model has its origin in compressive work on the matter 
of the core by gravitational forces. The trapped leptons 
are compressed and the matter is heated to high tem- 
peratures and thermal energies, both of which represent 
stored energy eventually to be radiated. Rotationally- 
powered and magnetic models ultimately derive their en- 
ergy from the conversion of gravitational binding energy 
changes during implosion into rotational kinetic energy 
(roughly conserving angular momentum), and then into 
magnetic energy. 

One of the characteristics of core-collapse supernovae 
that distinguishes them from thermonuclear (Type la) 
supernovae is that they leave a residue, the neutron star 
or black hole. It is not necessary to disassemble and un- 
bound this remnant to infinity, thereby paying a severe 
gravitational binding energy penalty. Since neutrino ra- 
diation renders the PNS more and more bound with time, 
if explosion were to require complete disassembly to in- 
finity, time would not be on the side of explosion. If fact, 
core-collapse supernovae would probably not be possible. 
However, in a fundamental sense, core-collapse supernova 
involve the transfer of energy from the core to the mantle, 
leaving the core behind. It is the mantle that is ejected. 
This mantle may start near ^100 km, not the canonical 
radius of ^10 km, and hence is much less bound. All 
CCSN explosion models are different models for core- 
mantle energy transfer, be it direct hydrodynamic (core 
piston), neutrino (mantle heating by core neutrinos), or 
MHD (mantle B-field amplification, tapping core rota- 
tion). One can bury in the residue a binding energy 
problem that could have gotten progressively worse with 
time. 

It may seem curious that the average thermonuclear 

supernova has an explosion energy that is similar to that 
of the average CCSN. However, the energy for the former 
derives from the burning of a large fraction of something 
like a Chandrasekhar mass, while an energy bound for the 
latter might be set in part by the gravitational binding 
energy of the stellar mantle surrounding a Chandrasekhar 
mass residue. The core-collapse explosion must eject this 
bound mantle. Burning yields ~0.5 MeV per baryon, and 
the binding energy of a Chandrasekhar mass is roughly 
nieC^ per electron, where nie is the electron mass. The 
latter obtains due to the fact that the Chandrasekhar 
mass is defined by the onset of relativity for the majority 
of its electrons. Since Ye is ^0.5, and perhaps 50% of 
a C/0 white dwarf burns to make a Type la supernova. 



both total energies are very approximately the number of 
baryons in a Chandrasekhar mass times ~0.5x0.5 MeV. 
This is very approximately one Bethe. The binding en- 
ergy of the white dwarf core of a massive progenitor and 
the binding energy of the stellar envelope around it will 
crudely scale with one another, due to the pseudo-power- 
law density profile of the latter. Therefore, the energy 
scales for both thermonuclear and gravitational super- 
novae explosions (the latter in the sense of a bound) are 
of comparable magnitude. This argument may be good 
to a factor of a few, but that it is good at all in a Universe 
with a much wider potential range of energies is perhaps 
noteworthy. 



1. Eigenvalue Problems 

There are two important approximate eigenvalue prob- 
lems associated with core collapse. The first involves 
the post-bounce, pre-explosion PNS structure, bounded 
by an accretion shock. The hydrodynamics during this 
phase is roughly quasistatic, so one can drop the time 
derivatives to arrive at a set of simultaneous ordinary 
differential equations for the hydrodynamic profiles inte- 
rior to the accretion shock. With shock outer boundary 
conditions, an inner core mass (Mc), a given accretion 
rate (M), and given core neutrino luminosities (Lj), one 
can convert this into an eigenvalue problem for the ra- 
dius of the shock (i?s)- One derives Rs in terms of the 
control parameters M, Li, and Mc, given assumptions 
about mantle heating due to Li and cooling due to elec- 
tron and positron capture on nucleons. As Burrows & 
Goshy (1993) showed, there is a critical curve in Lj ver- 
sus M space (ceteris paribus) above which there are no 
solutions to this eigenvalue problem. As one increases Li 
for a given M, Rg increases, but it can not increase to 
arbitrary values. At a critical curve Li, Rs is finite, but 
above the critical value of Lj for a given M the steady- 
state problem does not have a solution. The absence of 
a solution can be considered an approximate condition 
for explosion by the neutrino heating mechanism. The 
subsequent evolution is dynamical, with continued neu- 
trino heating depositing energy to power the explosion 
and expansion of the gas lowering the temperature (not 
the entropy!) and, thereby, the cooling rates. Further- 
more, as the mantle accelerates into explosion, the matter 
recombines from nucleons and alpha particles into iron 
peak nuclei (with {Z,A) depending upon Yg), thereby 
"returning" to the expanding matter the nuclear binding 
energy of the ejecta "lost" to shock photodissociation. 
These "original" ejecta may not contain more than a few 
hundredths of a solar mass, but since ^8—9 MeV are lib- 
erated per baryon, only ~0.1 solar masses of ejecta are 
needed to supply approximately one Bethe to the explo- 
sion. Similar critical curves can be derived that include 
multi-dimensional turbulence in 2D (Murphy & Burrows 
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2008) and 3D (Nordhaus et al. 2010) and these seem to 
be lower, facilitating explosion. Figure [TT] depicts the de- 
bris field of such a 3D explosion. However, the detailed 
reasons for this dimensional boost are still being studied. 

The second eigenvalue problem presents itself after ex- 
plosion and is the neutrino-driven wind that emerges 
from the proto-neutron star. An old model (Parker 1958) 
for the solar wind started by assuming that the plasma 
above the solar photosphere was in hydrostatic equilib- 
rium and that the energy luminosity due to electron con- 
duction through this atmosphere was constant. Since 
electron conductivity in a hot plasma depends almost 
solely on temperature (oc T^/^), one derives temperature 
as a function of radius (T cx l/r^/^). Hydrostatic equi- 
librium of an ideal gas can then be integrated to derive 
the pressure as a function of radius. What one finds is 
that around a spherical star the pressure must be finite 
at infinity! This shows that in order to maintain hydro- 
static equilibrium the star must be artifically embedded 
in a high pressure gas. Since the pressure in the inter- 
stellar medium is very low, this atmosphere cannot be 
stable and it would spontaneously erupt as a wind. The 
flow would transform into a steady-state outflow, with a 
sonic point and a supersonic asymptotic speed. Though 
we now know this particular physical model does not ap- 
ply to the Sun, these arguments led to the prediction 
of the existence of the solar wind. The wind mass loss 
rate would be a function of the driving luminosity, the 
stellar mass, and details of the heating process. It can 
be shown that if the derived temperature profile falls off 
more slowly than r~^, such an atmosphere is similarly 
unstable. 

The relevance of this scientific parable is that the same 
arguments can (more legitimately) be applied to the man- 
tle of the PNS. The balance of neutrino heating (oc l/r^) 
and neutrino cooling (cx T^) yields T oc l/r^/^, with 
a power-law index less than one. Therefore, without a 
bounding pressure the PNS atmosphere is unstable to a 
neutrino-driven wind (Duncan, Shapiro, & Wasserman 
1986; Burrows 1987; Burrows, Hayes, & Fryxell 1995). 
In the core-collapse context this bounding pressure is 
provided by the accretion ram and, while the shock is 
stalled, the wind is thwarted. However, after explosion 
and after the pressure around the PNS has subsided due 
to the progress of the supernova explosion, a neutrino- 
driven wind naturally emerges. The eigenvalue problem 
for M as a function of driving luminosity and PNS mass 
is easily solved. 

Therefore, in the context of the delayed neutrino heat- 
ing model, the supernova itself is the dynamical tran- 
sitional state between two quasi-steady-state eigenvalue 
problems, one accretion and the other a wind. The me- 
chanical power in the wind is lower than the instanta- 
neous power being poured into the early supernova be- 
cause the absorbing mass and neutrino optical depth (t^) 
of the atmosphere above the PNS around the gain region 



are much larger than in a tenuous wind. At the onset 
of explosion, how much larger the absorbing mass and 
depth are determines how much power is available for ex- 
plosion. If Ti, is large when the v^. and luminosities are 
large, their product will be large and the explosion will be 
robust. The simulated explosions of the 8.8-M0 progeni- 
tor model of Nomoto & Hashimoto (1988), with the very 
steep density gradient outside the core, transitioned so 
quickly into a wind that these model supernovae were ef- 
fectively wind-powered and had the correspondingly low 
explosion energy alluded to earlier (Kitaura et al. 2006; 
Burrows et al. 2007c). 



2. Simultaneous Accretion and Explosion 

When breaking spherical symmetry in the context of 
the multi-dimensional instabilities seen in modern super- 
nova simulations, a feature (some would say a virtue) 
of many proposed core-collapse supernova mechanisms 
is that during the early phases of explosion there can 
be simultaneous explosion and accretion (Burrows et al. 
2007a). Continued accretion onto the PNS from one di- 
rection can supplement the energy available to power ex- 
plosion in another. 

The neutrino mechanism, in part powered by accretion 
luminosity, is a good example. In spherical symmetry, ex- 
plosion is the "opposite" of accretion, and that source of 
neutrino driving subsides early after the onset of explo- 
sion. However, if the symmetry is broken and a neutrino- 
driven explosion first occurs in one direction, continued 
accretion onto the PNS from another direction can help 
maintain the driving neutrino luminosities. Though such 
accretion might be restricted to a small quadrant, neu- 
trino emissions are always much more isotropic than mat- 
ter distributions (Ott et al. 2008), with the result that 
accretion almost anywhere on a PNS surface leads to 
emitted neutrinos almost everywhere. In a sense, the 
same is true for the MHD mechanism, wherein explosion 
is bipolar along the rotation axis, while the spinning PNS 
accretes along the equator. Conserving angular momen- 
tum, such accreta continue to bring in the kinetic energy 
of differential motion needed to maintain the magnetic 
energy and pressure that power the explosion. As long 
as equatorial accretion continues, the core is an "engine" 
with a power source. After equatorial accretion ceases 
and the explosion assumes a more isotropic distribution, 
the engine subsides, but the supernova (or hypernova) is 
launched. The acoustic mechanism is the quintessential 
process that exploits simultaneous accretion on one side, 
which maintains the driving core g-mode oscillations, to 
power an explosion in the other direction (Burrows et al. 
2006,2007b). 
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3. Energetics 

Determining the energy of a detailed numerical explo- 
sion can be more awkward than one might think. Usually 
limited by the small size of the computational domain 
(e.g., 5000—20000 km), a successful shock wave encoun- 
ters this border and perforce stops within hundreds of 
milliseconds of the start of explosion and long before the 
explosion energy has asymptoted. For such calculations, 
neutrino energy deposition is still ongoing, recombina- 
tion of the nucleons and alphas has not completed, and, 
importantly, the baryon mass cut between the final PNS 
and the ejecta is not determined. In fact, the mass cut 
has never been consistently determined for any detailed 
numerical core-collapse supernova model. In addition, 
the explosion must work to unbind the star exterior to the 
mass cut, and this matter (most of the remaining star) 
can be bound by a few Bethes (see Figure |5|. The lar ger 
this binding energy, the stronger the explosion in the core 
needs to be to achieve a given final ejecta supernova ki- 
netic energy. The large magnitude of this binding energy 
for the most massive progenitors may be instrumental in 
either aborting what may have started as a promising su- 
pernova, or in ensuring that a black hole, rather than a 
neutron star, remains. In fact, the large binding energies 
for matter exterior to a canonical 1.5 Mq in progenitors 
that one has in the past thought should supernova and 
leave neutron stars (such as 20 or 25 Mq ZAMS stars) 
suggest that either the core explosion must be very vig- 
orous or that explosions in such stars fizzle. This would 
be unfortunate, since it is thought that the more massive 
progenitors are the likely primary sources for the oxy- 
gen that we see in abundance in the Universe. However, 
the relevant outer binding energies are functions not only 
of progenitor mass, but of modeler. This is yet another 
indication of the centrality of progenitor models to our 
understanding of the outcome of collapse. 

In the neutrino heating model, one way to achieve 
higher explosion energies may be to explode early. In the 
current paradigm, during the delay to explosion the neu- 
trino energy deposited is reradiated and useless. There 
is no accumulation of energy in the post-shock mantle 
until the explosion is underway. This suggests that a 
long delay to explosion may be detrimental, wasting as 
it does the neutrinos radiated by the core before the ex- 
plosion commences. However, an early explosion may be 
correspondingly useful, with perhaps some later fallback. 
Such an early onset may be easier in 3D (Nordhaus et al. 
2010). In addition, and counter- intuitively, both before 
and at the onset of explosion, the enthalpy fluxes and 
PNS mantle energies are negative, the latter at times 
even when the recombination energy is accounted for. 
For all the viable explosion mechanisms, the supernova 
does not attain its final energy at the instant of explosive 
instability, but must be driven after it starts and still 
needs to overcome the PNS and outer stellar envelope 



gravitational binding energies. Currently, no detailed 
neutrino-driving simulation has come within an order of 
magnitude of achieving this requirement, except perhaps 
the singular 8.8-Mq neutrino-wind-driven model. 



4. Conditions for Explosion by the Neutrino IVlechanism 

The neutrino mechanism, legitimately the front-runner 
in the CCSN mechanism sweepstakes, has engendered 
much speculation concerning the physical conditions for 
explosion. In section [IV. A. 1[ I described the critical con- 
dition between M, Li, and Mc that signals instability to 
explosion. This condition, with corrections to account for 
multi-D effects, still seems close to capturing the essence 
of the explosion condition. A detailed perturbation anal- 
ysis of these steady states to explore the complex eigen- 
frequencies of the monopolar and low-order pulsational 
modes, in particular to determine when their imaginary 
parts change sign, would lend a useful additional per- 
spective (Yamasaki & Yamada 2007). 

However, other explosion conditions have been prof- 
fered which aid understanding (e.g., Pejcha & Thompson 
2012; Burrows 1987; Janka 2001). By and large, all these 
are roughly equivalent. All sensible conditions must rec- 
ognize that since the matter interior to the pre-explosion 
shock is in sonic contact the condition for explosion must 
be a global one. A local condition, say on the pressure 
at the shock, has little meaning and can be misleading. 
It is the entire mantle structure that is exploding. More- 
over, the discussion in sectionllV.A.SIindicates how subtle 



things might be, with explosion commencing even when 
various otherwise obvious quantities associated with en- 
ergy or energy flux are negative. 

It should be noted that the gain region interior to the 
shock, in which neutrino heating outpaces capture cool- 
ing, surrounds a net cooling region where cooling dom- 
inates. The inner boundary of the cooling region coin- 
cides, more or less, with the and De neutrinospheres. 
This region gradually sinks in due to energy and lepton 
loss, steadily sending out rarefaction waves that are par- 
tially responsible for undermining the gradual outward 
progress of the quasi-static shock wave and the growth 
of the gain region. The larger the gain region and the 
smaller the cooling region the more likely the mantle is 
to explode. However, the cooling power generally out- 
strips net heating and this fact is one of the primary im- 
pediments to explosion. If there were no cooling region, 
or if the cooling power in the cooling region were signif- 
icantly reduced, neutrino absorption would quite easily 
lead to explosion. A focussed study on how Nature might 
accomplish this might bear fruit. 

A useful approximate condition for explosion is when 
the characteristic neutrino heating time in the gain re- 
gion (th) exceeds the advection time (Tadv) through it 
(Thompson et al. 2005; Janka 2001). For every set of 
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definitions of tfiese particular times, and there are a va- 
riety of definitions wliicli can vary by factors (Murpliy 
& Burrows 2008), the critical ratio itself should be cal- 
ibrated using hydro. In any case, r^dv can be set equal 
to Ar/veff, where Veff is some effective speed through 
the gain region that incorporates the sinuous trajectories 
of Lagrangian particles in multi-D. Quite naturally, Tadv 
is larger in 3D than in 2D, and larger in 2D than in ID 
(spherical). Tadv can also be written as AM/M, where 
M is the accretion rate through the shock and AM is the 
mass in the gain region. 

th can be defined as the internal energy in the gain 
region divided by the neutrino heating power. The latter 
is approximately L^^r^, where, again, Tj/ is the electron 
neutrino optical depth. Therefore, setting th equal to 
Tadv gives us Li,^ ^ {e/Ty)M, where e is the specific 
energy in the gain region, e might scale with the es- 
cape speed squared (We^c) ^-t the shock, and this quantity 
scales with the core mass and the inverse of the shock ra- 
dius (i.e., Ugg^ ~ 2GMc/Rs)- This rough relation yields a 
critical L,y^ versus M curve with a slope of (e/r^), which 
itself may be a weak function of L^^, M, and Mc that is 
better calculated numerically. However, this relation can 
be recast by noting that ^ n^^pAr, where k^^ is the 
electron neutrino absorption opacity and p is some mean 
mass density in the gain region. The result is 



(1) 



The actual constant of proportionality will depend upon 
details. Nevertheless, equation ([T]) states that the higher 
the absorptive opacity or the ratio of the mass in the 
gain region to the core mass the lower the critical lu- 
minosity for a given M and Rs- However, the quantity 
Rs/AM varies slowly with the control parameters L^^, 
M, and making equation([l]) a more direct connec- 
tion between them that succinctly summarizes the criti- 
cal condition of Burrows & Goshy (1993). A correction 
factor to account for multi-D effects could be added in 
the denominator. Not surprisingly, the critical curve and 
the Th ~ Tadv condition are roughly equivalent. 



5. Instability to Finite Perturbation 

The stalled accretion shock becomes unstable to out- 
ward expansion and explosion when (or near when) the 
critical L^^-M-Mc condition is met and exceeded. How- 
ever, it is also unstable to an abrupt, finite jump in its po- 
sition. If, by some mechanism, the shock were to be jolted 
suddenly to larger radii, the consequently lower matter 
temperatures would transiently result in a much lower 
integrated neutrino cooling rate behind the shock, while 
the larger radius of the shock would result in a larger gain 
region. Such a sudden, favorable, and finite shift to more 
heating and less cooling (and to a lower accretion ram 



pressure at the shock) could be irreversible and an ex- 
plosion might be ignited. However, the magnitude of the 
necessary finite perturbation is not known. The agency 
of such a jolt is also not known, but the accretion of den- 
sity discontinuities in the progenitor at composition and 
entropy boundaries (such as the inner edges of the silicon- 
or oxygen-burning shells) is seen in hydrodynamical sim- 
ulations to result in a quick outward (though modest) 
excursion in the average shock radius. If the actual den- 
sity jumps are larger than in the current generation of 
pre-coUapse models, or if there are significant variations 
in the density or velocity profiles of the post-bounce ac- 
creting matter, the necessary finite perturbations may 
be available. To be sure, this discussion is highly spec- 
ulative, but the possibility remains intriguing that the 
mantle shell interior to the stalled shock could be non- 
linearly unstable to explosion by large, finite-amplitude 
perturbations. 



B. Persistent Myths 

There are a number of what I would call "myths" that 
have emerged and persisted, despite compelling physical 
counter- arguments. To be sure, my list is idiosyncratic 
and the list of others may be different. One myth is that 
neutrino transport is more difficult than photon trans- 
port, requiring only a specialist's touch. In fact, non- 
LTE photon transport, with its multitude of level pop- 
ulations, spectroscopic data, coUisional processes, and 
lines, is much more difficult and challenging than neu- 
trino transport. The latter involves only continuum opac- 
ities and emissivities and one rate equation, that for Ye, 
which is coupled only to and i7g transport. True, there 
are six neutrino species and one does require knowledge of 
neutrino-matter couplings. Also, some astronomers shy 
away from the nuclear and neutrino realm, and are more 
interested in the dominant signatures of the Universe 
— those in photons. However, the physics of neutrino 
transport is far more straightforward than the physics 
of atomic and molecular spectroscopy and of the myriad 
coUisional processes in a heterogeneous soup of elements 
and ions. In addition, the numerical art of photon trans- 
port is a rigorous, well-developed, subject with expertise 
spread around the world, whereas experts in neutrino 
transport are few and far between. 

Another myth is that since ^-^3 x 10^'^ ergs of binding 
energy is emitted during the long-term ('-^ 10-50 second) 
PNS cooling and deleptonization phase, whereas the av- 
erage core-collapse supernova involves only ~10^^ ergs, 
the CCSN is a less than 1% affair, requiring exquisite 
precision and numerical care in approaching its theory. 
In fact, since the and 17^, absorption optical depths 
in the gain region are ~4— 10% and this is the fraction 
of the core Ve and i7e luminosities absorbed in the man- 
tle, the core-collapse neutrino mechanism is more like a 
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~4— 10% affair. Nevertfieless, it is often suggested that 
every detail makes a difference to tlie overall outcome, 
as if the mechanism itself hinged upon them. The result 
has been that minor effects have sometimes been allowed 
to loom large, often confusing those not intimating in- 
volved in the research. Examples are v-D annihilation, 
neutrino-electron scattering, electron capture on infall, 
neutrino-neutrino oscillations, and the nuclear symme- 
try energy, to name only a few. This is not to say that 
all these topics are not to be addressed, nor that the ul- 
timate theory can afford to ignore them. Incorporating 
the correct physics and performing detailed simulations 
will certainly be necessary to obtain the correct numbers. 
However, when the best, most detailed, extant exploding 
2D CCSN simulations may not be reproducing observed 
supernova energies by an order of magnitude perhaps a 
focus on details at the expense of global understanding 
is unfounded. Something much more important may be 
missing. 

A more innocuous myth is that a stellar-mass black 
hole can form directly. In fact, in the context of the col- 
lapse of an effective-Chandrasekhar-mass core, the inner 
homologous core will always rebound into the outer core, 
generating a shock wave. Interior to this shock wave at 
its inception and during its early life is only ~1.2— 1.5 
Mq of material and this is not enough for the core col- 
lectively to experience the general-relativistic instability 
that leads to stellar-mass black holes. Importantly, the 
inner core and shocked mantle are out of sonic contact 
with the supersonically infalling outer core, and, hence, 
do not yet "know" whether enough mass will accumulate 
to transition to a black hole. Sufficient matter must ac- 
crete through the shock before the core can go unstable. 
The wait, during which the core will fatten, might require 
hundreds of milliseconds to seconds, after which a second 
dynamical collapse to a black hole will ensue. Therefore, 
black hole formation is always preceded by an interme- 
diate PNS stage and cannot proceed directly (Burrows 
1984; Sumiyoshi, Yamada, & Suzuki 2007ab; Fischer et 
at. 2009; O'Connor & Ott 2011). 

Neutrino heating in the gain region naturally gener- 
ates a negative entropy gradient in the steady state. This 
gradient is unstable to convective overturn (the general- 
ized Rayleigh- Taylor instability) and leads to turbulence. 
This turbulence and the corresponding aspherical flow 
patterns, with neutrino-driven upflowing plumes and as- 
sociated downflows, has been seen in multi-D simulations 
since the 1990's and are central features of supernova the- 
ory. 

Later, Blondin, Mezzacappa, & DeMarino (2003) iden- 
tified an instability of the stalled shock, the standing ac- 
cretion shock instability ("SASI"), that has been well- 
diagnosed and studied by Foglizzo (2002,2009), Blondin 
& Mezzacappa (2006), Foglizzo, Scheck, & Janka (2006), 
Foglizzo et al. (2007), Yamasaki & Foglizzo (2008), 
Scheck et al. (2008), Iwakami et al. (2008), Kotake et al. 



(2009), and Foglizzo et al. (2012). In axisymmetric (2D) 
studies of the pure SASI (without neutrinos) , many wit- 
nessed a vigorous dipolar (^ = 1) "sloshing" mode that 
superficially seems like the corresponding sloshing motion 
seen in full 2D radiation/hydro simulations. This has led 
many, I believe incorrectly, to associate the motion seen 
in most full neutrino-transport runs with the motions 
seen in the simplified neutrino-free studies. Moreover, 
since the energy cascade in 2D turbulence (neutrino- 
driven or otherwise) is "backwards" (Boffetta & Musac- 
chio 2012), from small to large scales, and the SASI in- 
stability, by its nature, is on large scales (small spherical 
harmonic index £), this chance correspondence of domi- 
nant scales in 2D may also in part explain the confusion. 
More importantly, preliminary calculations performed in 
3D (e.g.. Burrows, Dolence, & Murphy 2012; Dolence, 
Burrows, & Murphy 2012; Hanke et al. 2012) do not 
show the "sloshing" dipolar motion along an axis many 
have come to associate with the SASI and that some have 
suggested is crucial to the CCSN mechanism (e.g., Marek 
& Janka 2009; Hanke et al. 2012). What are seen are 
bubble structures and plumes indicative of buoyant con- 
vection (Figure 12 1. Hence, the prominent £= \ "SASI" 



mode may be an artifact of 2D and its inverse cascade. 
However, when 3D simulations with reasonable physics 
become more readily available, the fact that the cascade 
in 3D is in the opposite direction, but the dominant SASI 
modes are still on large scales, should help clarify the true 
nature of the turbulence seen. 

I include this discussion on the SASI under "myths" 
because 1) current 3D simulations do not show the sig- 
nature features of the SASI ostensibly seen in 2D, and 
2) I fully expect that when credible and self-consistent 
3D radiation/hydrodynamic simulations are performed 
an objective analysis of its role will reveal it often to be 
subordinate to neutrino-heated buoyant bubble convec- 
tion. The reader is cautioned that not every supernova 
researcher will agree with this characterization, so caveat 
lector. Nevertheless, and curiously, though simulations 
with neutrino transport almost always show the classic 
rising bubble and downfiow patterns of buoyancy-driven 
non-linear convection, many people, even practitioners, 
had started to refer to all turbulent motions behind the 
shock as the "SASI." This may have served to confuse 
both insiders and outsiders alike. In summary, I am led 
to suggest that if neutrino driving is the energetic agency 
of explosion, it is naturally also the primary agency of the 
turbulence that aids explosion. 



V. CONCLUSIONS 

I have sought in this short review of core-collapse su- 
pernova theory to clarify its major facets, lay bear its 
physical underpinnings, and briefly summarize its current 
status, as I see it. Furthermore, I have tried to identify 
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some of what I consider to be myths that have crept into 
theoretical discourse. However, this colloquium is a per- 
sonal view of the theoretical landscape. There is clearly 
much more to be done before a cogent explanation of this 
central astrophysical phenomenon is available and veri- 
fied, and some of what I have suggested here may not 
survive future scrutiny. Nevertheless, it is hoped that at 
the very least this paper provides some novel and useful 
insights into the physics and astrophysics of both core 
collapse and supernova explosions. 
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FIGURES 
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FIG. 1 The electron fraction, Ye, versus enclosed mass for a suite of massive-star progenitor models from Woosley & Weaver 
(1995) ("sXXs") and Nomoto & Hashimoto (1988 — "heliumX" models). The numbers given (the "X"s in this caption) denote 
either the ZAMS masses or the helium-core masses. The dots indicate discontinuities and the three-pointed markers indicate 
the boundaries of fossil burning shells. The outer edges of the iron cores are well-marked by the large circular colored dots. 
Note that they are positioned from ~1.25 M© to ~2.0 Mq and are roughly in order of progenitor mass. See text for discussion. 
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FIG. 2 Similar to Figure [ij but depicting the logarithm of the mass density (p, in gm cm"'') versus interior mass (in solar 
masses) for various initial progenitor masses. Note that the lower-mass progenitors have steeper slopes, and that these profiles 
translate into more quickly dropping mass accretion rates (M) through the stalled shock after bounce. 
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FIG. 3 The velocity (in l(f cm s"'^) versus interior mass (in solar masses) at various pre- and post-bounce times. The 
shock wave, when present, is clearly indicated by the vertical drop, and is seen only for the post-bounce times. Within ^^20 
milliseconds of bounce, the shock wave has stalled into accretion and post-shock speeds go negative. Reviving this structure 
is the goal of modern core-collapse supernova theory. Note the line without the shock that changes slope in the middle. This 
transition in slope near ~0.6 Mq marks the edge of the homologous core a few milliseconds before bounce. Exterior to this 
minimum is the supersonic mantle whose maximum infall speed can reach ~80,000 km s~^. The 11-Mq progenitor model from 
Wooslcy & Weaver (1995) was used for this Figure, as well as for Figures 4, 6, 8, & 9. (Taken from numerical data generated 
in Thompson, Burrows, & Pinto 2003.) 
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FIG. 4 The luminosity (in units of 10^^ ergs s~^) at infinity of the i^e, Ve-, and, collectively, the z/^, Vr^ and Vr neutrinos, 
versus time (in seconds) around bounce. The Vc, breakout burst is clearly seen. After shock breakout, the temperatures are 
sufficient to generate the other species in quantity. Generally, the non-electron types carry away ~50% of the total, with the 
i^eS and z7eS sharing the rest. (Plot taken from Thompson, Burrows, & Pinto 2003.) 
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FIG. 5 The logarithm base ten of the gravitational binding energy (in ergs, including the thermal energy) of the shells in 
various progenitor stars (see Figure [T]) exterior to the interior mass coordinate (in M©) shown on the abscissa. These are the 
approximate energies that the supernova blast must overcome to eject the stellar shells exterior to a given residual neutron star 
or black hole mass. Note that the more massive the progenitor the greater the binding energy cost for a given baryon mass left 
behind. 



22 




FIG. 6 Snapshots of Ye versus interior mass (in solar masses) profiles before and just after bounce. Capture decreases Ye on 
infall, but lepton number is soon trapped in the interior. After bounce, the outward progress of the shock wave liberates 
neutrinos, creating a trough in Ye interior to the shock, but exterior to the opaque core. The shock wave resides just rightward 
of the steep drop in Ye on the right. The numerical data for this 11 Mq progenitor model run were taken from Thompson, 
Burrows, & Pinto (2003). 




FIG. 7 A storyboard of the evolution of the core of the "onion-skin" structure into the radiating proto-neutron star. See text 
for a discussion. 
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FIG. 8 Temperature (in MeV) versus interior mass (in solar masses) at various times before and just after bounce. The shock 

is depicted by the vertical drop in temperature and propagates out in mass as (and after) it stalls into accretion. Initial 
post-bounce central temperatures arc ~10 MeV and the peak temperature just before shock breakout can exceed 20 MeV, but 
soon falls. The various neutrinophcrcs after ~20— 100 milliseconds after bounce are at ~1. 0—1.1 M©. (Numerical data taken 
from Thompson, Burrows, & Pinto (2003).) 
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FIG. 9 Similar to Figure|8] but for entropy (in ks per baryon) versus interior mass (in solar masses). The black curves depict 
the evolution of profiles using realistic neutrino transport and stopping ~I0 milliseconds after bounce, while the red curves 
depict the corresponding developments when neutrino physics is turned off. Note that the red curves extend further in mass 
and reach much higher early mantle entropies. See text for discussion. (Numerical data taken from Thompson, Burrows, & 
Pinto (2003).) 
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FIG. 10 An early snapshot during the bipolar explosion of the mantle of a rapidly-rotating core. Depicted are representative 
magnetic field lines. The scale is 2000 km from top to bottom. The extremely twisted lines are just interior to the shock wave 
being driven out by magnetic pressure generated within 200 milliseconds of bounce in a "2. 5" -dimensional magneto-radiation- 
hydrodynamic calculation conducted by Burrows et al. (2007d). The progenitor employed for this calculation was the IS-M© 
model of Heger et al. (2005). 
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FIG. 11 The debris field generated in a 3D neutrino-driven explosion approximately 200 milliseconds after its onset. The scale 
from top to bottom is 1000 km. The blue exterior is a rendering of the shock wave, the colored interior is a volume-rendering of 
the entropy of the ejecta, and the sphere in the center is the newly-born neutron star. (Numerical data taken from Nordhaus 
et al. (2010) and the IS-M© progenitor from Woosley & Weaver 1995 was used.) 
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FIG. 12 A difTerent rendering of the 3D neutrino-driven explosion shown in Figure [TT] but approximately 25 miUiseconds after 
its onset. The scale from top to bottom is ~500 km. The prominent bubble structures in magenta are isoentropy surfaces and 
are provided to highlight the neutrino-heated bubbles that seem to be generic in current 3D simulations. The bounding shock 
is not shown. The red colored interior (between the bubbles) is an experimental (only partially successful) volume-rendering 
of the density. (Numerical data taken from Nordhaus et al. (2010) and the 15-M0 progenitor from Woosley & Weaver 1995 
was used.) 



